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ABSTRACT 

We present a new SCUBA image of the cluster MS 0451.6— 0305, which exhibits strong, 
extended sub-mm flux at 850 /im. The most striking feature in the map is an elongated region 
of bright sub-mm emission, with a flux density of ~ lOmJy over several beam-sizes. This 
region is apparently coincident with a previously known optical arc (which turns out to be 
a strongly lensed Lyman Break Galaxy at z = 2.911), as well as with a newly identified 
multiply imaged ERO (Extremely Red Object) pair predicted to be at a similar, if not identical 
redshift. By combing a detailed lensing model with deep images from HST, Chandra, CFHT, 
JCMT, and spectra from the VLT, we conclude that both the strongly lensed optical arc and 
ERO systems have properties consistent with known sub-mm emitters. Using a simple model 
for the two sources, we estimate that the multiply lensed EROs contribute the majority of the 
flux in the SCUBA lensed arc. Correcting for the lensing amplification, we estimate that the 
inherent 850 fluxes for both objects are < 0.4 mJy. If the LBG and ERO pair are truly at 
the same redshift, then they are separated by only ~ lOkpc in the source plane, and hence 
constitute an interacting system at z ~ 2.9. Higher angular resolution observations in sub- 
mm/mm will permit us to more accurately separate the contribution from each candidate, and 
better understand the nature of this system. 
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1 INTRODUCTION 



Since the installation of the Submillimetre Common User Bolome- 
ter Array (SCUBA; Holland et al. 1999) at the JCMT telescope, 
hundreds of luminous dusty galaxies have been detected (see Blain 
et al. 2002 for a review). Such observations are still far from rou- 
tine, however, with the detection rate in random fields being typi- 
cally one or two sources of modest signal-to-noise ratio (SNR) per 
shift of telescope time. Bright sources, which are easier to obtain 
detailed sub-mm/mm follow-up observations for using heterodyne 
receivers or bolometers tuned to other wavelengths, are correspond- 
ingly rarer. This, coupled with the fact that the dust responsible 
for the sub-mm luminosity absorbs radiation at other wavelengths, 
means that SCUBA detected galaxies are often extremely faint in 
the optical. Hence, redshifts, morphologies and spectral energy dis- 
tributions (SEDs) have proven elusive for most sources. 

One technique which can help, pioneered by Smail, Ivison 
& Blain (1997), is to use the lensing amplification of rich galaxy 
clusters to boost the detectability of background sub-mm galax- 
ies. The rate at which sources are detected is increased by a fac- 
tor ~ 3, with an improved prospect of discovering the occasional 
bright object at relatively high SNR. Here we present detailed ob- 
servations of just such a source in the field of the massive cluster 
MS 045 1.6—0305, which we discovered in an earlier survey of 9 
rich cluster fields (Chapman et al. 2002a, hereafter C02). This pa- 



per uses new SCUBA data, as well as rigorous chop deconvolution, 
in order to produce a map free from the chopping artifacts common 
in sub-mm observations. 

The sub-mm emission in this image is extended along a giant 
optical arc discovered in this cluster by Luppino et al. (1999, here- 
after L99), and thus there is little doubt that the SCUBA emission 
is lensed as well. In this paper, we combine deep optical, near-IR 
(NIR), and X-ray data with a detailed lensing model in order to 
interpret the sub-mm emission in this remarkable SCUBA field. 



2 OBSERVATIONS 

MS 045 1.6—0305 is the most luminous X-ray cluster in the Ein- 
stein Observatory Extended Medium Sensitivity Survey (EMSS; 
Gioia et al. 1990). It is a 2 = 0.55 cluster, which has been stud- 
ied extensively in the optical (Gioia & Luppino 1994, hereafter 
L99) and X-Ray (Molnar et al. 2002, Donahue et al. 2003; here- 
after D03). It makes an ideal SCUBA target because it is relatively 
compact compared to the SCUBA array size, and the strong radio 
Sunyaev-Zerdovich (SZ) decrement (Reese et al. 2000) suggests 
that the SZ increment at 850 fim might be detectable. The SZ and 
X-ray signals clearly indicate that this is a massive cluster, and the 
presence of a giant optically detected arc unambiguously verifies 
that it is lensing background galaxies. Thus SCUBA observations 
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have the potential to detect a combination of diffuse SZ emission, 
background lensed sources, and galaxies in the cluster itself. 



2.1 Sub-mm observations 

SCUBA 'jiggle-mode' maps of the region were taken simultane- 
ously at 850 /im and 450 fim in two separate JCMT runs; The first 
was on 03 September 1998 and the second was 10-11 Novem- 
ber 1999. Each has approximately the same observing time of 
~ 22 ksec, but the weather in 1999 was significantly better, with an 
average atmospheric opacity of rsso — 0.25, compared to 0.40 in 
1998. In addition, upgrades after the 1998 run improved SCUBA's 
overall sensitivity by roughly 10% at 850 fim and 50% at 450 fim. 

To remove the strong atmospheric emission component of the 
data, the JCMT secondary mirror 'chops' between the target and 
two nearby reference positions on the sky. Since the sky emission 
is common to all positions, but the astronomical target is not, the 
difference between the data at these positions removes the major- 
ity of the common mode noise. Three different chop directions, all 
using 30 arcsecond throws, were used in 1998. Analysis of these 
data, revealed that the region was excessively rich in bright sub-mm 
sources, raising concern that the chopping may have removed some 
of the measured signal. The sub-mm map of MS 0451.6—0305 pre- 
sented in C02 use these data exclusively. 

In an effort to detect the Sunyaev-Zel'dovich (SZ) increment 
in the cluster, the 1999 observations were taken with very wide 
chop throws (180 arcsecond), so that signal within the 2 arcminute 
SCUBA field of view would not be chopped out. Note that chop- 
ping onto sources outside the field is unlikely to be a concern, since 
the chop was held at constant azimuth and therefore sky rotation 
would quickly move sources out of the chopped position. Zemcov 
et al. 2003 use the 1999 data only, since the data in C02 used as 
a chop throw too small to be useful for SZ studies. By comparing 
the maps made from each set, it is clear that 1998 data has, in fact, 
chopped out some signal inside the field. 

The map presented here combines these two datasets, and 
given the differences between the 1998 and 1999 maps, it seems ap- 
propriate to be fairly comprehensive in describing the map-making 
procedure which we adopted. This 'direct-inversion' method is de- 
scribed in the appendix. It could be adapted for all SCUBA jiggle- 
maps which have in-field chopping. However, it is much more com- 
putationally intensive than standard approaches, and will only lead 
to mild improvements for general 'blank' fields. On the other hand, 
for fields with multiple or extended sources, a method of this sort 
is essential in order to recover the correct fluxes across the map. 

Because the 1999 data are roughly three times more sensitive, 
the gains in co-adding the 1998 data are modest. Nevertheless the 
maps from the 2 separate runs are consistent with each other af- 
ter the deconvolution is performed. We present our best 850 fim 
and 450 fim maps in Fig. Q The 850 fim image shows a continu- 
ous ridge of sub-mm emission, which was partly removed in the 
chopped data presented in C02. More importantly, the act of chop- 
ping in-field shifted the centroid south-westward by removing flux 
on the edge of the arc region. The 450 fim image, on the other hand, 
contains no significantly detected sources. The lack of a detection 
of the primary arc with the 450 fim array (3er upper limit of 48 mJy) 
places a lower limit of z ~ 2 for the source redshift, for reasonable 
dusty galaxy SEDs. We also note that since the 1999 data were 
taken with a 16-point jiggle pattern (thus fully sampling at 850 fim, 
but not at 450 fim), the resulting 450 fim map has quite inhomoge- 
neous noise. 



2.2 HST Observations 

MS 045 1 .6-0305 was imaged with the F702W filter on WFPC2 in 
1995 as part of a program to study luminous X-ray clusters (project 
ID 5987, PI: Donahue). We down-loaded the processed, calibrated 
image directly from the archive at the Canadian Astronomy Data 
Centre. More recently, the Advanced Camera for Surveys (ACS; 
Ford et al. 1998) was used to obtain deep images with the F775W 
and F850LP filters (project ID 9292, PI: Ford). We obtained the 
raw data from the public archive and reduced them using the MUL- 
TIDRIZZLE package (Koekemoer et al. 2002). 

The WFPC2 F702W images are composed of 4 independent 
dithered exposures, with a total exposure time of 10.4 ksec. In each 
of the two ACS bands, two exposures totaling 2.4 ksec were ob- 
tained. However, because these observations were not dithered, the 
resulting ACS images are prone to detector artifacts and cosmic ray 
contamination. This makes estimating fluxes for some objects in the 
field problematic. Astrometry for the WFPC and ACS images was 
fixed using 13 objects from the GSC2.2 catalogue. The 5a limiting 
magnitudes for the F702W, F775W, and F850LP images are 26.3, 
24.9 and 25.1 AB magnitudes, respectively. 

The HST F702W image is presented in Fig. [2] and shows in 
detail the lensed arcs reported originally in L99. ARC1 appears 
to contain knots surrounded by diffuse emission which were not 
previously resolved from ground-based imaging. These knots dis- 
play a mirror symmetry, arguing that ARC1 corresponds to a fold- 
arc (two images merging together). We use this fact later in creat- 
ing a lensing model for the cluster. North-east of the cluster cen- 
tre is a second arc: ARC2, showing a very smooth light distribu- 
tion. VLT spectra described later show that ARC 1/2 are at very 
different redshifts, and thus are not two images of the same ob- 
ject. Furthermore, the colours of these two arcs are quite different 
(mF702\v — J71F850LP = 0.3 AB for ARC1 versus 0.6 for ARC2), 
a finding first reported in L99 using lower resolution V- and R- 
band images. 

Another notable object in the field is a very blue 
(m,F702w — rriFg,s,OLP = 0.0 AB), morphologically disturbed 
source (which we label 'P') that seems to be coincident with 
the SCUBA point source eastward of the arc. We note also that 
the Chandra contours extend around this source, suggesting the 
presence of another massive component to the cluster. 



2.3 Near IR Observations and EROs 

The CFHT-IR camera on the 3.6-m Canada-France-Hawaii Tele- 
scope (CFHT) was used to image MS 0451.6—0305 in each of the 
J-, H-, and A"'-band filters. The images were obtained in moder- 
ate conditions, with seeing roughly 1.0 arcsec. In each filter, ob- 
servations were dithered in a hexagonal pattern, providing uniform 
coverage of a 3 x 3 arcminute region. Several short (< 90 second) 
integrations were taken, for a total exposure time of 7.2 ksec in 
H, 5.6 ksec in K' , and 6.8 ksec in J. Observations of the standard 
FS11 were taken to measure the zero points. The limiting magni- 
tudes in the J-, H- and if'-band images were 24.4, 23.5, and 23.8 
(AB, 5(j) respectively. 

Data were reduced using the DIMSUM package in IRAF, to- 
gether with SEXTRACTOR to detect objects and estimate fluxes 
within 2 arcsecond apertures. We compared the _R"'-selected cat- 
alogue against 84 objects detected in the shallower survey of Stan- 
ford et al. (1998). Although there is no evidence of bias between 
our catalogues, we did find an RMS scatter of 0.1 magnitudes be- 
tween sources in each of the three NIR bands, and we use this as 
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Figure 1. SCUBA maps of MS 0451.6—0305. The 850 (im flux contours at 4, 6, 8, and lOmJy are overlaid on each. The approximate FWHM of the beam-size 
at each wavelength is shown as the circle in the bottom right corner of each map. The 850 /im image clearly shows emission extended over several beams, 
while the short wavelength image shows no discernible flux. 




Figure 2. HST F702W image of MS 0451.6— 0305. LEFT: This deep, high resolution image clearly shows the 2 gravitationally lensed arcs previously 
identified in lower resolution images (L99). The stronger arc to the east (Arc 1) appears to be made up us a series of 'knots'. Its counter-image, ARC1 
ci' is also circled. Overlaid are the SCUBA 850 (ira flux contours with the same levels as in Fig.0 The dotted contours are taken from the medium-band 
(1.5-4.5 keV) Chandra image, with levels at 0.03, 0.05, 0.1, 0.2 and 0.3 photons cm" 2 pixel -1 (pixel size is 0.492" in the Chandra image). Positions of 
some of the red objects detected in the NIR images are also overlaid using circles. X-ray point sources from Molnar et al. (2002) coincide with two of these 
objects, and are shown as crosses. We also circle the blue point source 'P', which is coincident with the eastern sub-mm point source. TOP RIGHT: A close 
up of ARC1 from the F702W image. BOTTOM RIGHT: A close up of the ARC1 counter-image (denoted by the arrow). The circle denotes the position of 
B3/C3. Both close-up images are roughly 4" wide. 
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the lower limit to the error estimate for our J-, H-, and K'-band 
magnitudes. Our magnitudes also agree with a deeper Subaru NIR 
image described in Tanaka et al. (2003). 

A colour-magnitude diagram (Fig.[5) using the J and K' mag- 
nitudes, and a false colour image of the cluster created using the 
NIR images (Fig. |4j reveals objects with extremely red colours 
(J — K' > 2.3). In Section[3]we will describe the lensing model 
of the cluster, but it is important here to highlight sources B and 
C, which are a pair of EROs re-imaged in 3 places due to strong 
lensing. The first and second image, sources Bl/Cl and B2/C2 are 
mirrored across the critical line in the lens model, and have simi- 
lar magnification factors. Unfortunately B2/C2 are in a particularly 
crowded region of the image, making accurate photometry difficult. 
At the position of B3/C3, the lensing is much weaker, and thus the 
object separation decreases to the point where the two distinct ob- 
jects overlap. 

The connection between sub-mm galaxies and some EROs is 
well established; dust responsible for the sub-mm emission also 
reddens the optical-UV spectrum. NIR bright objects invisible in 
deep HST imaging were found coincident with sub-mm sources 
in observations towards other lensing clusters (Smail et al. 1999). 
Roughly one-third of the sub-mm galaxies detected in the UK 8- 
mjy survey have an ERO as the most likely counterpart (Ivison 
et al. 2002). In addition, Webb et al. (2004) claim that pairs of 
EROs denote regions of sub-mm emission, though not necessar- 
ily the source of it. And in general, stacking analyses of optical and 
NIR selected galaxies in sub-mm fields show a convincing trend for 
redder objects to be statistically more SCUBA-bright (Wehner et al. 
2003, Borys et al 2004). These findings, along with the fact that the 
3 ERO pairs here trace out the sub-mm contours quite convincingly, 
make this a strong candidate for the sub-mm arc counterpart. 

2.4 Chandra observations 

MS 045 1.6-0305 was observed with the Chandra Advanced CCD 
Imaging Spectrometer (ACIS) detector for a total of 41ksec. Full 
processing details are contained in D03. The astrometry for a Chan- 
dra ACIS-S observation is typically accurate to about 1.0", but can 
have errors as large as 3.0". X-ray point source fluxes are described 
in Molnar et al. (2002). Two sources are detected within the field 
of the SCUBA observations, and are coincident with CI and C3. 
The X-ray source coincident with CI was detected with 45.2 ± 8.1 
counts, while that of C3 has 21.4 ± 5.7. The ratio of counts is the 
same as the ratio of magnification factors derived from the lensing 
model described in the next section. 

Since a point source was not reported in D03 at the position 
of C2, we estimated the X-ray flux using a 2" aperture. The lensing 
model predicts that that CI and C2 have approximately the same 
amplification factor, thus we expect ~ 45 counts. However, we 
only find 13 ± 6 counts between 0.5-2.0 keV at the position of C2. 
We note that photons from C2 pass through more of the cluster gas, 
since it is closer to the brighter X-ray centre than CI (see Fig.|2j- 
If we assume that the NIR fluxes are accurate, then the fainter C2 
(by 0.5 magnitudes in K') suggests the number of counts from C2 
should be around 25. Our estimate of 13 counts is still low com- 
pared to this expectation, but given the difficulty in extracting an 
X-ray flux estimate of point sources on top of a brighter, diffuse 
background (the cluster gas), it does not seem unreasonable. More 
accurate NIR imagery, coupled with more detailed modeling of ex- 
tinction through the cluster, could resolve this issue. Based on the 
strength of the lensing model and colours/geometry of the EROs, 
we still conclude that the three B/C pairs are the same objects, with 
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Figure 3. Colour-magnitude diagram of MS 0451.6— 0305 using the NIR 
data. The region above the dashed line is beyond the limit of the survey. 
Only objects within 1' of the BCG are plotted. Labeled sources denote the 
triply-imaged ERO pair described in the text. None of the points have been 
corrected for lensing amplification. All fluxes were measured in 2" aper- 
tures, except for B2, which required a 1" aperture due to its proximity to 
nearby sources in the image. 



C also being an X-ray source (despite the apparent non-detection 
of X-rays in C2). 

The X-ray profile discussed in D03 allows us to predict the SZ 
increment at 850 /im, which, in the absence of chopping, peaks at 
roughly 3 mjy. By applying the astrometry from the sub-mm data, 
we constructed a simulated map of what the SZ signal would look 
like when the chopping is accounted for. The peak drops to 2 mjy 
and the SZ emission falls off to roughly 1 mjy at the position of 
the main arc. If we subtract this shape from the data it makes little 
difference to the flux levels or SNRs in the SCUBA map. On the 
other hand, the contamination of so much lensed sub-mm emission 
makes this particular field useless for sub-mm SZ studies, at least 
until much higher resolution images can be obtained (Zemcov et 
al. 2003). Indeed, measuring the SZ increment in any strong lensing 
cluster will prove difficult simply because the resolution is too poor 
to remove the enhanced number of point sources detected in the 
field. 



2.5 VLT spectrum of ARC1 

Deep multi-object spectroscopy of MS 045 1.6—0305 was obtained 
in service mode for the FORS2-GTO programme: 68.A-0015 (PI: 
Appenzeller) on December 10th, 11th, 12th and 21st 2001, with 
the FORS-2 instrument (Appenzeller & Rupprecht 1992) on the 
VLT-Yepun (UT4). We retrieved these data from the ESO/VLT 
archive and reduced them using standard IRAF procedures. Two 
masks were prepared ('M128' and 'M129') and exposed with the 
high throughput/low resolution grism GRIS_150I (5.2A/pixel). For 
each mask five exposures of 2.9 ksec each were obtained (totaling 
14.5 ksec). ARC1 is part of the 'M129' mask, and ARC2 of the 
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Figure 4. A true-colour image of MS 0451.6—0305 using the NIR data. The scale is approximately the same as in Fig.|2] The three images of the ERO pair 
are labeled here, as well as the blue source 'P' (although it is essentially invisible in the NIR). The outer white contour traces shows the z = 2.911 critical 
tangential line, and the inner one corresponds to the radial line. Sub-mm contours at the same levels in Fig. 2 are shown as the dashed lines. 



'Ml 28' mask. The four nights had variable seeing, ranging from 
0.6 to 1.0", and were likely not fully photometric. Nevertheless we 
obtained a crude flux calibration of these observations using the 
spectrophotometric standard star GD108. 

ARC2 displays only one strong emission line at A = 7144A. 
We identify it as Oil, and thus calculate a redshift of z — 0.917 ± 
0.001. The derived redshift for ARC1 (Fig.|3J is 2.911 ± 0.003 
based on the strongest interstellar absorption lines: Sill A1260.42A, 
OIA1302.17A and CIIA1334.5A as well as CIVA1548.19 and 
1550.77 in absorption. There is also a hint of the triple SiVI lines at 
~ A1397A and the Sill A1526.70A line. 

The ARC1 spectrum lacks AGN signatures sometimes found 
in the spectra of radio selected sub-mm galaxies. Moreover, it ap- 
pears to show Ly q in absorption, a feature seen in only about 20% 
of the spectroscopic survey of ~ 100 sub-mm galaxies conducted 
by Chapman et al. (2004). That ARC1 does not look similar to 
most known sub-mm galaxy spectra does not necessarily mean it is 
not the host of the sub-mm emission. The set of blank-field sub- 
mm galaxies with spectroscopically determined redshifts are bi- 
ased towards galaxies with strong Ly a emission, and it is clear we 
are studying a sub-mm galaxy which is intrinsically much fainter 
(and possibly much different) than those reported in Chapman et 
al. (2003, 2004). 



3 A LENSING MODEL FOR MS 0451.6 0305 

In our first attempt to model the system we employed a simple 
isothermal ellipsoid for the cluster mass (using parameters from 
D03), and a second elliptical mass distribution to model the Bright- 
est Cluster Galaxy (BCG). This allowed us to gauge the general 



properties of the lens, and showed that the ridge of sub-mm emis- 
sion generically lies along the critical line for high redshift sources. 
However, the importance of determining the precise position and 
magnifications of the multiply lensed EROs necessitated a more de- 
tailed approach. Indeed, as the Chandra contours in Fig.|2|show, the 
mass within the central 1' of the cluster is not smoothly distributed, 
and the BCG position and X-ray centroid are not coincident. 

To model the mass distribution of this cluster we used both a 
cluster mass scale component and cluster galaxy mass component 
in a similar way to Kneib et al. (1996; see also Smith et al. 2003). 
Using the LENSTOOL software (Kneib 1993) we interactively im- 
plemented the lensing constraints. The high resolution HST data 
show a clear mirror symmetry along the ARC1, making it the merg- 
ing of 2 images (fold arc). A simple elliptical model allows us to 
identify its counter-image 3" to the south-east of B3/C3, which is 
much less amplified (see Table0. 

Similarly, the ERO B/C pair can then be included in the lens- 
ing model. Bl/Cl and B2/C2 are straddling the critical line at red- 
shift z — 2.9, thus by placing the B/C pair at z ~ 2.9, we easily 
identify B3/C3 as the third image of the B/C pair. 

By combining the ARC1 and B/C constraints, we refined 
the mass model, and are able to place more precise limits of 
z = 2.85 ± 0.10 for the ERO B/C pair. It is therefore likely that 
both ARC1 and B/C are at the same redshift. If it is the case, then 
the separation between the LBG and the ERO pair in the source 
plane is only ~ lOkpc, suggesting that these 3 galaxies are inter- 
acting. 

Using the photometry derived from the optical and near- 
infrared data, as well as the magnifications estimated from the lens- 
ing model, we summarize, in TableQ the list of objects highlighted 
in Fig.|2|and Fig.|4] 
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Figure 5. FORS2 spectrum of ARC1 plotted in the observers rest frame. The bottom curve is the raw spectrum, and a smoothed version is plotted above it. 
Based on the strong absorption features, the redshift of the arc is constrained to be z = 2.911 ± 0.003. Ly a absorption is evident near 4766A, which is seen 
in only about 20% of the submillimetre galaxy spectra in the sample of Chapman et al. (2004). 



Table 1. Photometry for objects of interest in the field of MS 0451.6—0305. These limits are shown as 5(7, and all are reported in the AB magnitude system. 
None of the fluxes are corrected for lensing magnification. All magnitudes are measured in 2" apertures, except for the arcs (measured using ellipses covering 
the same area in each of the fields, with the counter-image ARClci measured using a 1" aperture because of possible contamination from a nearby source) and 
B2, where a magnitudes in a 1" aperture is also given because of possible flux contamination. The final column gives magnifications from the lensing model; 
Most are well constrained, except for ARC1. Since it lies right along the critical line, an more accurate estimate of the amplification is difficult. 



ID 


RA 


DEC 


F702W 


F775W 


F850LP 


J 




H 




K< 




J-K' 


Mag. 


ARC1 


04:54:13.0 


-03:00:55.9 


22.9 ±0.1 


22.6 ±0.1 


22.6 ±0.1 


21.2 ±0 


1 


21.0 ±0 


1 


20.2 ±0 


1 


1.0 


5CM00 


ci(l") 


04:54:11.1 


-03:01:26.6 


25.7 ±0.4 


25.2 ±0.2 


25.3 ± 0.3 


> 24.4 




> 23.5 




> 23.8 






4.7±0.5 


Bl 


04:54:12.8 


-03:00:44.1 


24.9 ±0.3 


25.0 ±0.2 


24.4 ±0.2 


23.1 ±0 


1 


21.4 ±0 


1 


20.6 ±0 


1 


2.5 


8±1 


B2 


04:54:12.7 


-03:01:01.3 


> 26.3 


> 24.9 


> 25.1 


22.9 ±0 


1 


22.2 ±0 


1 


20.7 ±0 


1 


2.2 


10±1 


B2(l") 






> 26.3 


> 24.9 


> 25.1 


23.8 ±0 


2 


22.9 ±0 


2 


21.3 ±0 


1 


2.5 




CI 


04:54:12.8 


-03:00:47.6 


> 26.3 


> 24.9 


> 25.1 


> 24.4 




> 23.5 




21.6 ±0 


1 


> 2.8 


10±1 


C2 


04:54:12.7 


-03:00:59.1 


> 26.3 


> 24.9 


> 25.1 


> 24.4 




> 23.5 




22.1 ±0 


1 


> 2.3 


11±1 


B3/C3 


04:54:10.9 


-03:01:24.5 


25.2 ±0.2 


> 24.9 


> 25.1 


24.1 ±0 


2 


21.8 ±0 


1 


21.0 ±0 


1 


3.1 


5±1 


P 


04:54:15.6 


-03:00:38.7 


24.0 ±0.2 


23.9 ±0.1 


24.0 ±0.2 


> 24.4 




> 23.5 




> 23.8 






1.5 



4 DISCUSSION 

With the lensing model in hand, we are now able to further exam- 
ine which galaxy hosts the sub-mm emission. Apart from the (we 
believe remote) possibility that the sub-mm emission arises from 
other sources entirely, there are 3 obvious possibilities, which we 
discuss in turn. We also comment on the additional sub-mm source 
to the east of the arc, which we label 'P'. 



4.1 Scenario 1: The blue arc as the source of the sub-mm 
emission. 

Although the spectrum for ARC1 does not match that of a 'typi- 
cal' sub-mm galaxy from the Chapman et al. (2004) sample, it is 
similar to the most absorbed quartile of the Lyman break galax- 
ies (LBGs), which generally correlate with the reddest, most dust 
extincted varieties (see Shapley et al. 2003). It is this dust which 
can be responsible for the sub-mm emission, though using the rest 
frame UV spectra of LBGs to estimate the sub-mm flux is notori- 
ously difficult (Chapman et al. 2000). In a sample of ~ 30 Lyman 
break galaxies with the most extreme reddening (and hence most 



likely chance to be sub-mm bright), only one (MMD-11) was de- 
tected. 

Indeed, several other studies find a similar lack of sub- 
millimeter emission from individual LBGs (Webb et al. 2003, Bo- 
rys et al. 2004). On the other hand, Peacock et al. (2001) stack 
LBGs in the deep sub-mm map of the Hubble Deep Field (Hughes 
et al. 1998) to estimate that the average 850 fim flux from an LBG 
is 0.20 ± 0.04mJy/M Q /i _2 yr _1 . Given that the unmagnified sub- 
mm flux would be 0.1 — 0.2 mjy if the LBG is the host, then it 
is entirely possible that we have simply detected an ordinary LBG 
producing stars at a rate of a few solar masses per year. 

An argument against ARC1 being the sole sub-mm source is 
the fact that the sub-mm emission is much more extended than the 
optically detected arc. It extends at least 15 arcseconds further north 
(a full SCUBA beam-width). Over 45 arcseconds (3 beams) sepa- 
rate the arc from B3/C3, where the toe of the sub-mm emission to 
the south-west tapers off. Note that the ARC1 counter-image is too 
faint to be responsible for the toe, since the amplification factors 
predicted from the lens model are not consistent with the observed 
sub-mm flux. Even assuming the lowest reasonable magnification 
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factor for ARC 1 , we would predict only ~ 1 mjy of 850 fim flux at 
the position of the ARC1 counter-image. 



4.2 Scenario 2: Contribution of the multiply lensed EROs to 
the sub -mm flux. 

Galaxies with very red colours often suggest the presence of large 
amounts of dust, and several sub-mm blank field surveys (Ivison et 
al. 2002, Webb et al. 2003, Borys et al. 2004) have found such ob- 
jects coincident with sub-mm detections. The triply-imaged ERO 
pair in MS 045 1.6—0305 is therefore a reasonable candidate for 
the counterpart to the sub-mm arc. The relative fluxes of objects 
Bl/Cl and B3/C3 match the relative amplifications predicted from 
the lensing model. B2 is a little difficult to compare, given the un- 
certainty in the flux due to the contamination from a nearby source, 
and the discrepancy with C2 has already been discussed. 

Unlike the last scenario, the position of the three images can 
easily accommodate the extremities of the extended emission. The 
sub-mm contours around B3/C3 suggest that it is a SCUBA point 
source, though the signal-to-noise is insufficient to be certain. As- 
suming B/C is the host of the sub-mm emission, and using the peak 
850 fim flux of 12mJy, we would predict 6mJy at the position of 
B3/C3. This is essentially what is observed. If the EROs are the host 
of the sub-mm emission, then the unmagnified flux is ~ 1.0 mjy at 
850 ^m. 

It is also worth noting the several red objects in southern part 
of field. These are in the same region as the southernmost extension 
of the sub-mm arc, and might be contributing sub-mm flux: Using 
colour-colour plots and the Pozzetti & Mannucci (2000) criteria for 
distinguishing dusty starbursts from passively evolved ones, Tanaka 
et al. (2003) show that both object B and an red object 4" to the 
northeast of B3/C3 are consistent with dusty starbursts. If so, this 
would help explain the extended sub-mm emission between B3/C3 
and its counter-images. 

In Fig.|5|we plot the optical fluxes for the best detected object 
from this trio of pairs, Bl, on a template SED obtained from the 
library of Efstathiou, Rowan-Robinson, & Siebenmorgen (2000). 
These combine the Bruzual & Chariot stellar population synthesis 
models (Bruzual & Chariot 1993) with a radiative transfer code , 
and produce SED templates that match local star-burst galaxy spec- 
tra. We translated the templates to z — 2.85 and scaled them to 
match the NIR fluxes. As Fig.|6|shows, there is reasonable agree- 
ment between the assumed redshift, the measured photometry, and 
the SED template. There is some freedom in these models, and 
hence we have not attempted to provide detailed constraints on the 
star-burst age, formation history, and composition. The main point 
is that the data support the claim that the EROs can be responsible 
for the sub-mm flux. Ground based NIR spectroscopy, as well as 
Spitzer observations, will be of considerable help in determining 
the nature of this galaxy pair. 



4.3 Scenario 3: Contribution from both 

We have now established that both objects are capable of producing 
sub-mm flux, so naturally we need to ask if both objects are con- 
tributing to the sub-mm arc. To test this, we created a blank map the 
same size as the SCUBA data, and inserted sources at positions cor- 
responding to Bl, B2 and B3, as well as ARC1. The relative fluxes 
between B 1, B2 and B3 was fixed by the lensing model predictions 
for amplification, leaving the flux ratio between ARC 1 and B 1 as 
the only remaining parameter. The final constraint is that the peak 




Observed frame wavelength/yttm 



Figure 6. Spectral Energy Distribution deduced from the HST and NIR 
data for the ERO Bl (circles) and ARC1 (triangles). The square point is 
the 850 fim emission of sub-mm arc. The 5a upper limits is also shown for 
the 450 fim measurement. The solid line is a model template redshifted to 
z = 2.85, and scaled to match the NIR measurements of Bl. The dotted 
line is another model redshifted to z = 2.911 and scaled to match the 
NIR data from ARC1. These models are derived from Efstathiou, Rowan- 
Robinson, & Siebenmorgen (2000), and assume a 72Myr old star-burst in 
the case of Bl, and 45 Myr for ARC1. Though the NIR fluxes are close fits 
to the model by design, the measured 850 /im flux is only a factor of 2 from 
the predictions. 

flux needs to match that observed, which in this case is 12 mly. The 
results of these simulations are shown in Fig.0 

We found that no contribution from ARC1 resulted in B 1/2 be- 
ing resolvable, which does not match the observed sub-mm image. 
In addition, ARC 1 needs to contribute at least 2 mly in order to 'fill 
in' the gap between B 1/2. At more than 7 mly, ARC 1 dominates the 
contours, and fails to reproduce the extended nature of the source. 
Although it is difficult to precisely constrain the relative contribu- 
tions, this exercises suggests that sub-mm emission from both the 
LBG arc and the ERO pair are required to explain the morphology 
and brightness of the sub-mm arc, with ~ 2/3 of the observed flux 
coming from the ERO B/C pair. 



4.4 The sub-mm point source, 'P' 

Though the ERO pair and LBG arc just discussed are the focus of 
this paper, it is worth noting the blue galaxy 'P', which seems coin- 
cident with the SCUBA detected point source north-east of the arc. 
Because there is simply no other optically detected galaxy within 
the error circle (6") of the SCUBA source, it is reasonable to con- 
clude that the objects are at least associated with each other, even 
if they are not necessarily the same thing. 

Note in Fig.|2|how the X-ray contours bend around this object. 
This suggests the presence of another massive component interact- 
ing with the cluster, or perhaps a projection effect from another 
X-ray cluster. If the sub-mm object is behind the cluster, then it 
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Figure 7. Simulating the relative contributions of the LBG and ERO. The 
plus symbol in each panel denotes the centre of ARC1, and the axes are in 
units of arcseconds. In the top left panel we show the 850 (im contours at 4, 
6, 8 and 10 mJy from the original sub-mm image. The top right shows a sim- 
ulation where a source was placed at the position of Bl, B2, and B3 (with 
relative fluxes determined from the lensing model), but no flux from ARC1. 
The bottom left shows a model where the LBG arc contributes 4 mJy, and 
then 8 mJy in the last panel. These simulations suggest that that at least 
some flux is necessary from both ARC 1 and B/C in order to reproduce the 
observed sub-mm image. 

is only weakly lensed; the prediction is ~ 1.5 from our lensing 
model, meaning that the unmagnified 850 /im flux for this object 
would be about 4 mJy. The lack of a 450 fim detection is not very 
constraining. Using a wide range of temperatures and emissivities 
in a single dust component SED, we predict a 450 pirn flux < 48 for 
z > 0.5. Since the 450 /im map has a 3a upper limit comparable 
to this, we cannot rule out that the sub-mm source is actually com- 
ing from the cluster itself. This would be an interesting case, since 
it would suggest that star-formation is occurring in a massive sub- 
component merging with MS 0451.6—0305. Further spectroscopic 
and X-ray studies may unravel this particular mystery. 



would place them only ~10 kpc apart at z ~ 2.9, and implies that 
the objects are interacting and their interaction is likely to be the 
origin of a violent star-burst revealed by the strong sub-mm emis- 
sion. To confirm that the LBG is effectively spatially connected 
with the ERO pair can only be determined by measuring the red- 
shift of the EROs using NIR spectroscopy or via millimeter CO 
line search; Given the constraints on the redshifts of either candi- 
date, sub-mm CO observations can be used to detect the CO(3-2) 
transition line that falls within the observation window of several 
ground-based heterodyne receivers. 

To verify these conclusions, and to accurately separate the 
contributions from each component will require sub-mm/mm in- 
terferometry. The SCUBA flux is bright enough that it should be 
detectable using the IRAM Plateau de Bure interferometer. The 
SMA and ultimately ALMA will also have the power to resolve 
the components in great detail. With such further study, and using 
the lensing amplification of nature's telescope, we should be able 
to provide a full accounting of the source of sub-mm emission in 
this intriguing system. 
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5 CONCLUSIONS AND FUTURE WORK 

We have presented broad-band sub-mm and near-IR images of the 
cluster MS 0451.6— 0305 and attempted to understand the nature 
of the highly extended sub-mm emission. The only other case of 
a multiply lensed sub-mm detection in a cluster field has recently 
been reported in Abell 2218 (Kneib et al. 2004). They identified 
the source of the triple sub-mm image as a LBG galaxy (showing a 
complex structure, indicative of a close merger) in a z ~ 2.5 group. 

In the case of MS 0451.6—0305, the evidence suggests that the 
sub-mm emission is coming from at least two components: a blue 
arc at z — 2.911 and a multiply imaged ERO pair predicted from 
the lens model to be at z ~ 2.85. Using a simple model for the 2 
sources, we estimate that the multiply lensed EROs contribute ~2/3 
of the observed flux in the sub-mm arc. Correcting for the lensing 
amplification, we estimate that the unlensed fluxes at 850 ^im are 
0. 1 mJy for the LBG and 0.4 mJy for the ERO pair. 

The uncertainty in the redshift for the ERO keeps the possibil- 
ity open that the ERO pair and LBG are at identical redshifts. This 
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APPENDIX A: MAP-MAKING BY DIRECT MATRIX 
INVERSION 

To make maps from data taken with SCUBA, one normally adds 
the measured signal into the pixel corresponding with the position 
of a bolometer in the 'on' position. For observations with in-field 
chopping however, this will result in a map with negative echos. 
Since the MS 0451.6—0305 field appears to be quite dense in bright 
sub-mm sources, it is necessary to deconvolve the chop to obtain a 
single-beam map. This can be posed as a linear algebra problem. 
Following the notation given in Stompor et al. (2002), the M mea- 
surements that make up the time-stream data vector, dj are related 
to a pixelized region of sky, mp that has elements via a sparse 
M x N mapping matrix, A: 



9 

d t = .4m p + n t (Al) 

The vector nt is the noise estimate for each measurement with sta- 
tistical properties described by M =< ntn t < >. The best estimate 
of the flux in the pixelized map, mp, in a least squares sense, has a 
well known solution (Stompor et al. 2002, Tegmark 1997, Wright 
1996): 

m p = (AN" 1 A 7 ')' 1 AATMt, (A2) 
M = (AAf- 1 A T y 1 . (A3) 
M is an N x N matrix that measures how much the pixels 

1/2 

are correlated with each other. M.J is the formal error estimate 
for pixel i. Ais typically extremely sparse, with each row here hav- 
ing only three non-zero entries [1.0, —0.5, —0.5] corresponding to 
the pixels that the primary and 2 off bolometers hit. Given that 
SCUBA observations are differential in nature, the matrix AN^A 
is formally singular, since there exists a non-zero vector m p where 
.4mp = 0. Specifically, and unsurprisingly, the degeneracy is in 
the overall mean of the map, which SCUBA is insensitive to (i.e., 
setting m p to a constant). 

J\f measures how correlated the noise is in the time domain. 
For anything but uncorrelated noise this matrix is difficult to man- 
age without special techniques, since N is a M x M matrix 
and M can easily reach 10 6 samples. The combined dataset for 
MS 0451.6-0305 has 1.14/0.52 x 10 6 samples at 450/850 ^m , 
corresponding to 12 hours of integration time. 

Assuming uncorrelated noise, the matrix AN~ X A and vector 
AJV~ 1 d[ can be populated in time order, requiring O(M) opera- 
tions. Thus we are faced only with the task of inverting an N x N 
matrix. 850 fim jiggle maps are typically binned into 3" pixels in 
order to reasonably sample the beam. Given the 2.3' field of view 
of a jiggle map observation, the number of pixels is on the order 
of a few thousand, which is small enough that any modern work- 
station can invert the required matrix and obtain the solution in an 
acceptably short time. 

The only caveat is in the selection of 'dummy pixels', which 
are used to account for positions a bolometer may hit that are out- 
side the solution region. For example, the 1999 data used 180 arc- 
sec chop-throws, which place the bolometers in the off position 
well outside the central target region. The procedure we adopt is 
to assume that there is no signal outside the target region, and 
assign a single dummy pixel. To verify the procedure, simulated 
data-sets were created, based on the astrometry information in the 
MS 0451.6— 0305 observations. False sources were placed such 
that some of the time the chop would cancel out their signals. To 
correct for the mean of the map (which is forced to zero), we select 
a region of the map free of sources and compute the average signal, 
which is then applied to the whole map (there is some subjectivity 
to this procedure, since for real maps one will not know for certain 
which regions are devoid of signal). The \ 2 f° r the reconstructed 
maps compared against the input source map was always within 
la of the number of pixels in the map. In doing these tests, we also 
tested varying the number and position of dummy pixels to verify 
that they did not adversely affect the inversion. 

This algorithm was then applied to the actual 
MS 0451.6—0305 data. First, maps from the 2 separate runs 
were solved and compared. The \ 2 between these half maps is bet- 
ter than that obtained from the raw maps. The overall improvement 
was not very dramatic but this is not completely unexpected since 
the chop configuration only causes 2 of the sources in the map 
to interfere with each-other. Also, the 1999 data are of superior 
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quality, and hence is weighted higher in the inversion. However, 
the improvement in consistency between the half maps is dramatic 
in the regions where there are chopped sources. 

Although the gain may be marginal in some cases, in general 
this direct inversion technique can be straightforwardly applied to 
SCUBA jiggle-map data.lt is easy to automate the procedure, with 
the only required input being how to treat bad bolometers, how to 
assign pixels outside the target region, and what to do about the DC 
level. 



